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Abstract 

A topological insulator protected by time-reversal symmetry is realized via spin- 
orbit interaction driven band inversion. The topological phase in the Bii_ x Sb x sys- 
tem is due to an odd number of band inversions. A related spin-orbit system, the 
Pbi ^Sn^Te, has long been known to contain an even number of inversions based on 
band theory. Here we experimentally investigate the possibility of a mirror symmetry 
protected topological crystalline insulator phase in the Pbi^^Sn^Te class of materi- 
als which has been theoretically predicted to exist in its end compound SnTe. Our 
experimental results show that at a finite-Pb composition above the topological in- 
version phase transition, the surface exhibits even number of spin-polarized Dirac 
cone states revealing mirror-protected topological order distinct from that observed 
in Bii x Sb x .. Our observation of the spin-polarized Dirac surface states in the inverted 
Pbi ^Sn^Te and their absence in the non-inverted compounds related via a topolog- 
ical phase transition provide the experimental groundwork for opening the research 
on novel topological order in quantum devices. 
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Topological insulators are band insulators featuring uniquely spin-polarized surface states, 
which is believed to be a consequence of odd number of inversions of bulk bands driven by 
spin-orbit interaction alone or in combination with crystal lattice modulations Q-Q. The 
odd number of band inversion transitions in the bulk demarcate a Z 2 topological insulator 
phase with spin-polarized surface states from that of a conventional band insulator phase 
as previously demonstrated in topological insulator Bi 1 _ 2 .Sb :r and the Bi 2 Se 3 series [l-5]. 
Within band theoretical calculations, lead tin telluride Pbi-^Sn^Te, which is a narrow band- 



aa, 



gap semiconductor widely used for infrared optoelectronic and thermoelectric devices 
has long been known to contain an even number of band inversions as electronic structure is 
tuned via the Sn/Pb ratio js 13] . Theoretical calculations have also predicted the occurrence 
of surface or interface states within the inverted band-gap upon the band inversion transition 



llNl3| but no experimental evidence of the surface states has been reported so far. Recent 
theoretical works [l^ . [lol ] have further stimulated the experimental search for surface states 
in Pbi-^Sn^Te. In a recent calculation fl5|. Hsieh et al. predicted that the topological 
surface states in SnTe |l5| are protected by crystalline space group symmetries, in contrast 
to the time-reversal symmetry protection in the well known Z 2 topological insulators. 

The crystal group symmetry of Pbi-^Sn^Te, critical for the realization of the topological 
crystalline insulator (TCI) phase Ijj], is based on the sodium chloride structure (space group 
Fm3m (225)). In this structure, each of the two atom types (Pb/Sn, or Te) forms a separate 
face-centered cubic lattice, with the two lattices interpenetrating so as to form a three- 
dimensional checkerboard pattern. The first Brillouin zone (BZ) of the crystal structure 
is a truncated octahedron with six square faces and eight hexagonal faces. The band-gap 
of Pbi_ x Snj;Te is found to be a direct gap located at the L points in the BZ [8]., which 
are the centers of the eight hexagonal faces of the BZ. Due to the inversion symmetry of 
the crystal, each L point and its diametrically opposite partner on the BZ are completely 
equivalent. Thus there are four distinct L point momenta. It is well established that the 
band inversion transitions in the Pbi-^Sn^Te take place at these four L points of the BZ 

9Hl3|. As a result, even number of (four) inversions exclude the system from being a Kane- 
Mele Z 2 topological insulator |2| . However, it is interesting to note that the momentum-space 
locations of the band inversions coincident with the momentum-space mirror plane within 
the BZ. This fact provides a clue that band inversions in Pbi-^Sn^Te may lead to a distinct 
topologically nontrivial phase that is irrelevant to the time-reversal symmetry but may be 
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the consequence of the spatial mirror symmetries of the crystal, which has been theoretically 



shown in the inverted end compound SnTe 



151 ] . However, it is experimentally known that 



SnTe crystals actually subject to a rhombohedral distortion [16], which, strictly speaking, 
breaks the crystal mirror symmetries and hence excludes a gapless TCI phase 15]. More 
importantly, SnTe crystals are typically heavily p-type [n| due to the fact that Sn vacancies 
are thermodynamically stable 17| . which makes the chemical potential of SnTe to lie below 
the bulk valence band maximum, consequently not cutting across the surface states. The 
surface states, which may exist within the band-gap, are thus unoccupied and cannot be 
experimentally studied by photoemission experiments 18]. On the other hand, the Pb-rich 
samples are reported to possess the ideal (averaged) sodium chloride structure without global 



rhombohedral distortion 
the predicted TCI phase 



gl which thus preserve the proper crystal symmetry required for 



15]. Moreover, it is possible to grow Pb-rich crystals, in which 



the chemical pot ential can be brought up above the bulk valence band maximum (in-gap or 
even n-type) 19] so as to lie near the surface state Dirac point. 

Therefore, in order to experimentally investigate the possibility of a topological phase 
transition, and thus rigorously proving the existence of a TCI phase in the Pbi-^Sn^Te 
system, we hereby utilize angle-resolved photoemission spectroscopy (ARPES) and spin- 
resolved ARPES to study the low energy electronic structure below and above the band 
inversion topological transition in the Pb-rich compositional range of the system, through 
which the observation of spin-polarized surface states in the inverted regime and their ab- 
sence in the non-inverted regime is demonstrated, correlating the observed spin-polarized 
Dirac surface states with the expected band inversion topological phase transition. We fur- 
ther map out the critically important spin structure of the surface states protected by mirror 
symmetries in the inverted topological composition. These results show that the topological 
order observed in Pbi_ x Sn x Te is distinct from that previously discovered in Bi!_ x Sb x or 
Bi 2 Se3 [l-^]. Our spin-resolved study of Pbi_ x Sn a .Te paves the way for further investigation 
of many novel properties of topological crystalline order 14], I2CN22I] in real materials and 
devices. 



Results 

Comparison of non-inverted and inverted compositions 

In order to capture the electronic structure both below and above the band inversion 
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transition (theoretically predicted to be around x ~ 1/3 0,[io|), we choose two representative 
compositions, namely x = 0.2 and x = 0.4, for detailed systematic studies. Figure [lb shows 
the momentum-integrated core level photoemission spectra for both compositions, where 
intensity peaks corresponding to tellurium 4d, tin 4d, and lead 5c? orbitals are observed. 
The energy splitting of the Pb orbital is observed to be larger than that of the Sn orbital, 
which is consistent with the stronger spin-orbit coupling of the heavier Pb nuclei. In addition, 
the spectral intensity contribution of the Sn peaks in the x = 0.4 sample (red) is found to 
be higher than that of in the x = 0.2 sample (blue), which is also consistent with the larger 
Sn concentration in the x = 0.4 samples. We perform systematic low energy electronic 
structure studies on these two representative compositions. Since the low energy physics 
of the system is dominated by the band inversion at L points (L points projected onto X 
points on the (001) surface), we present ARPES measurements with the momentum space 
window centered at the X point, which is the midpoint of the surface BZ edge (see Fig. [Tb). 
Figure Eb,b show the ARPES Fermi surface and dispersion mappings of the Pbo.8Sno.2Te 
sample (2 = 0.2). The system at x — 0.2 is observed to be gapped: No band is observed 
to cross the Fermi level in the Fermi surface maps (Fig. Eb). The dispersion measurements 
(Fig. [2b ) reveal a single hole-like band below the Fermi level. This single hole-like band 
is observed to show strong dependence with respect to the incident photon energy (see 
Fig. [2b, and also detailed in Supplementary Figures S1-S3 and Supplementary Methods 
in the supplementary information (SI).), which reflects its three-dimensionally dispersive 
bulk valence band origin. As a qualitative guide to the ARPES measurements on x = 0.2, 
we present first-principles based electronic structure calculation on the non-inverted end 
compound PbTe (Fig. |2b). Our calculations confirm that PbTe is a conventional band 
insulator, whose electronic structure can be described as a single hole-like bulk valence band 
in the vicinity of each X point, which is consistent with our ARPES results on Pb-rich 
Pbo.8Sno.2Te. The three-dimensional nature of the calculated bulk valence band is revealed 
by its k z evolution in Figure [2b, which is in qualitative agreement with the incident photon 
energy dependence of our ARPES measurements shown in Figure [2b. 

Now we present comparative ARPES measurements under identical experimental condi- 
tions and setups on the Pbo.6Sno.4Te (x = 0.4) sample. In contrast to the conventional band 
insulator (insulating) behavior in the x = 0.2 sample, the Fermi surface mapping (Fig. [2H) 
on the x = 0.4 sample shows two unconnected metallic Fermi pockets (dots) on the opposite 
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sides of the X point. Such two-pockets Fermi surface topology cannot be easily interpreted 
as the bulk valence band since the low energy bulk valence band of the Pbi-^Sn^Te sys- 
tem is a single hole-like band, unless certain surface-umklapp processes of the bulk states 
are considered, which usually have only weak cross-section in ARPES measurements. The 
assignment of the observed two-Fermi-pockets as surface-umklapp processes or other bulk- 
related origins can be further ruled out by incident photon energy dependence [k z dispersion) 
measurements and our spin polarization studies of these metallic states (see Supplementary 
Figures S1-S5 in the SI). The dispersion measurements on the x = 0.4 sample are shown in 
Figure [2fe. The single hole-like bulk valence band, which is similar to that in the x = 0.2 
sample, is also observed below the Fermi level. More importantly, a pair of metallic states 
crossing the Fermi level on the opposite sides of the X point is observed along the f — X — T 
mirror line momentum space direction. These states are found to show no observable dis- 
persion upon varying the incident photon energy (further details in Supplementary Figure 
S3), reflecting its two-dimensional character. On the other hand, the single hole-like band is 
observed to disperse strongly upon varying the incident photon energy, suggesting its three- 
dimensional character. At a set of different photon energy (k z ) values, the bulk valence 
band intensity overlaps (intermixing) with different parts of the surface states in energy 
and momentum space. At a photon energy of 18 eV, the intermixing (intensity overlap) is 
strong, and the inner two branches of the surface states are masked by the bulk intensity. 
At photon energies of lOeV and 24 eV, the surface states are found to be relatively better 
isolated. These ARPES measurements suggest that the x = 0.4 sample lie on the inverted 
composition regime and that the observed surface states are related to the band inversion 
transition in Pbi-^Sn^Te as predicted theoretically llMl3l. [15|]. As a qualitative guide, we 
present first-principles based electronic structure calculation on the inverted end compound 
SnTe (Fig. [2f). The calculated electronic structure of SnTe is found to be a superposition of 
two k z nondispersive metallic surface states and a single hole-like k z dispersive bulk valence 
band in the vicinity of the X point, which is in qualitative agreement with the ARPES 
results on Pbo.6Sno.4Te. 



Surface state topology of Pbo.6Sno.4Te 

We now perform systematic measurements on the surface electronic structure of the 
Pbo.6Sno.4Te. Figure [3b shows the wide range Fermi surface mapping covering the first 
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surface BZ. The surface states are observed to be present, and only present, along the mir- 
ror line (r — X — r) directions. No other states are found along any other momentum 
directions on the Fermi level. In close vicinity to each X point, a pair of surface states 
are observed along the mirror line direction. One lies inside the first surface BZ while the 
other is located outside. Therefore, in total four surface states are observed within the first 
surface BZ, in agreement with the fact that there are four band inversions in Pbi_ x Sn x Te. 
The mapping zoomed-in near the X point (Fig. [3b) reveals two unconnected small pockets 
(dots). The momentum space distance from the center of each pocket to the X point is 
about 0.09 A . Dispersion measurements (Ep vs k) are performed along three important 
momentum space cuts, namely cuts 1, 2, and 3 defined in Fig. [3b, in order to further reveal 
the electronic structure of the surface states. Metallic surface states crossing the Fermi level 
are observed in both cuts 1 and 2, whereas cut 3 is found to be fully gapped, which is consis- 
tent with the theoretically calculated surface states electronic structure shown in Figure [Tfe. 
In cut 1 (Fig. [3bl), which is the mirror line (F — X — T) direction, a pair of surface states 
are observed on the Fermi level. The surface states in our Pbo.6Sno.4Te samples are found to 
have a relatively broad spectrum, which can be possibly understood by the strong scattering 
in the disordered alloy system similar to the broad spectrum of the topological surface states 
in the Bi^^-Sb^ alloy 23] . In addition to the scattering broadening, the surface states are 
also observed to tail on the very strong main valence band emission (for example, the white 
region of intensity distribution in Fig. [3bl). In the case of Figure [2b (photon energy of 24 eV), 
the bulk valence band maximum (VBM) locates outside the first surface BZ. The surface 
state inside the first surface BZ is relatively better isolated from the bulk bands as compared 
to the one outside the first surface BZ. We thus study the dispersion along cut 2 (Fig. [3bl), 
which only cuts across the surface states inside the first surface BZ. Both the dispersion 
maps and the momentum distribution curves in cut 2 reveal that the surface states along 

ermi level. In many topological 



24j . Fitting of the momentum 



cut 2 are nearly Dirac-like (linearly dispersive) close to the 
insulators, the surface states deviate from ideal linearity 
distribution curves of cut 2 (see Supplementary Figure S6 for data analysis) shows that the 
experimental chemical potential (Ep) lies roughly at (or just below) the Dirac point energy 
(Ed), Ep = Ep, ± 0.02 eV. The surface states' velocity is obtained to be 2.8±0.1 eV-A 
(4.2 ±0.2 xlO 5 m/s) along cut 2, and l.ldz0.3 eV-A (1.7 ±0.4 xlO 5 m/s) for the two outer 
branches along cut 1, respectively. 
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Spin polarization measurements of the surface states 



We study the spin polarization of the low energy states of the Pbo.6Sno.4Te samples, 
which are highly dominated by the surface states near the Fermi level. We further compare 
and contrast their spin behavior with that of the states at high binding energies away from 
the Fermi level, which are highly dominated by the bulk valence band in the same sample. 
Spin-resolved (SR) measurements are performed in the spin-resolved momentum distribution 



curve mode 



25 



26], which measures the spin- resolved intensity and net spin polarization 



at a fixed binding energy along a certain momentum space cut direction (detailed in the 
Methods section). As shown in Figure Hh, our spin-resolved measurements are performed 
along the mirror line (T — X — T) direction, since the electronic an d sp in structure along 
this direction is most critically relevant to the predicted TCI phase [l5j. Considering that 
the natural Fermi level of our x = 0.4 samples are very close to the Dirac point (which 
is spin degenerate), the spin polarization of the surface states are measured at 60 meV 
below the Fermi level in order to gain proper contrast, namely SR-Cut 1 in Figure UK- As 
shown in the net spin polarization measurement of SR-Cut 1 in Figure 0b, in total four spins 
pointing in the (±) in-plane tangential direction are revealed for the surface states along 
the mirror line direction. This is well consistent with the observed two surface state cones 
(four branches in total) near an X point along the mirror line direction. To compare and 
contrast the spin polarization behavior of the surface states (SR-Cut 1) with that of the bulk 
states, we perform spin-resolved measurement SR-Cut 2 at Eb = 0.70 eV, where the bulk 
valence bands are prominently dominated. Indeed, in contrast to SR-Cut 1 reflecting the 
surface states' spin polarization, no significant net spin polarization is observed for SR-Cut2, 
which is expected for the bulk valence bands of the inversion symmetric (centrosymmetric) 
Pbx_3.Snj.Te system. Our observed spin polarization configuration of the surface states is also 
in qualitative agreement with the first-principles calculation spin texture of the SnTe TCI 
surface states (see Supplementary Figures S7-S9 for texture calculation). And Experimental 
derived topological invariant (Mirror Chern number 27J, |28j) tim = —2 also agrees with 



theoretical prediction for SnTe 



15|. 



Discussion 

We discuss the possibility that the observed surface states in our data are the signature of 
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the theoretically predicted TCI phase. The TCI phase was predicted to be observed in the 



band inverted side and absent in the non-inverted side [15j. In our data of the Pbo.6Sno.4Te 
samples which lie on the inverted side, surface states are observed, and only observed, along 
the two independent mirror line (F — X — F) directions. Within the first surface BZ, two 
surface states are observed on each mirror line, which are found to locate in vicinity of the 
X points. Dispersion measurements (e.g. Cut 2 in Fig. [3]) reveal a nearly linear dispersion. 
Spin polarization measurements show that these surface states are spin-polarized and their 
spin polarization direction is locked with their momentum (spin-momentum locking). The 
overall electronic structure observed for our surface states are in qualitative agreement with 
the theoretically predicted surface states of the TCI phase in SnTe 15[. It is important to 
note that in an alloy, the space group symmetry is only of the averaged type. The spatial 
mirror symmetry required for the topological protection is only preserved in a globally 
averaged sense. To this date, it is not even theoretically known whether only averaged mirror 
symmetry is sufficient for the predicted TCI phase. In the absence of any theoretical work for 
the alloy, the experimental proof of mirror protection in the alloy Pbi_ x Sn :r Te(Se) system 
perhaps requires the observation of strictly gapless surface states. The current ARPES 



works 



29 



30] have not demonstrated resolution better than 5 meV, and thus fitting of the 



energy-momentum distribution curves lacks the resolving power to exclude gap values less 
than 5 meV (demonstrated in Supplementary Figure S10 and Supplementary Discussion). 
On the other hand, a proof of strictly gapless surface states in fact requires a surface and 
spin sensitive experimental setup with energy resolution better than 1 meV, which can even 
exclude possible small gap values of < 1 meV at the Dirac point. Even a gap size of even less 
than 1 meV can lead to observable change in spin-transport experiments, thus destroying 
the delicate TCI phase. Therefore, while the existence of surface states correlated with 
band inversion transition is established here in our work, their gapped or gaplessness nature 
requires surface transport experiments on samples with improved quality, which is currently 
not possible 29]. 

Our observations of spin-momentum locked surface states in the inverted composition and 
their absence in the non-inverted composition suggest that Pbi-^Sn^Te system is a tunable 
spin-orbit insulator analogous to the BiT^Si^Se^ system j^], which features a topological 
phase transition across the band inversion transition. However, the key difference is that 
odd number of band inversions in the BiT^Si-aSe^ system leads to odd number of surface 
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states, whereas even number of band inversions in Pbi-^Sn^Te system leads to even number 
of surface states per surface BZ. We further present a comparison of the Pb .6Sn 4 Te and a 
single Dirac cone Z2 topological insulator (TI) GeBi2Te4 3lM33l|. As shown in Figure EK-c, 
for GeBi2Te4, a single surface Dirac cone is observed enclosing the time-reversal invariant 
(Kramers') momentum T in both ARPES and calculation results, demonstrating its Z2 
topological insulator state and the time-reversal symmetry protection of its single Dirac 
cone surface states. On the other hand, for the Pb .6Sn .4Te samples (Fig. EH-f), none of the 
surface states is observed to enclose any of the time-reversal invariant momentum, suggesting 
their irrelevance to the time-reversal symmetry related protection. With future ultra-high- 
resolution experimental studies to prove the strict gap less nature of the Pbo.6Sno.4Te surface 
states and therefore the topological protection by the crystalline mirror symmetries, it is 
then possible to realize magnetic yet topologically protected surface states in the Pbi-^Sn^Te 
system due to its irrelevance to the time-reversal symmetry related protection, which is 
fundamentally not possible in the Z 2 topological insulator systems. Considering the wide 
applications of magnetic materials in modern electronics, such topologically protected surface 
states compatible with magnetism will be of great interest in terms of integrating topological 
insulator materials into future electronic devices. 

We note that one advantage of the Pbi_ x Sn x Te system is that it can be easily doped with 



34-361. 



manganese, thallium, or indium to achieve bulk magnetic or superconducting states 
The symmetry in the Pbi^Sn^Te system (its nonmagnetic character) can be broken by 
magnetic or superconducting doping into the bulk or the surface. In future experiments, 
it would be interesting to explore the modification of our observed surface states brought 
out by magnetic and superconducting correlations, in order to search for exotic magnetic 
and superconducting order on the surface. Such magnetic and superconducting orders on 



the surface states in 
topological insulators 



bj T Sn T Te can be different from those recently observed in the Z 2 



37 



due to its very distinct topology of surface electronic structure. 
The novel magnetic and superconducting states to be realized with this novel topology 
are not strongly related to the question of gapless or gapped nature of the TCI phase. 
Therefore our observation of the spin-polarized surface states presented here provides the 
much desired platform for realizing unusual surface magnetic and superconducting states in 
future experiments. 
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Methods 

Electronic structure measurements. Spin-integrated angle-resolved photoemission 
spectroscopy (ARPES) measurements were performed with incident photon energies of 8 
eV to 30 eV at beamline 5-4 at the Stanford Synchrotron Radiation Lightsource (SSRL) in 
the Stanford Linear Accelerator Center, and with 28 eV to 90 eV photon energies at beam- 
lines 4.0.3, 10.0.1, and 12.0.1 at the Advance Light Source (ALS) in the Lawrence Berkeley 
National Laboratory (LBNL). Samples were cleaved in situ between 10 to 20 K at chamber 
pressure better than 5xl0 -11 torr at both the SSRL and the ALS, resulting in shiny surfaces. 
Energy resolution was better 15 meV than and 1% of the surface BZ. In order to cross-check 
the cleavage surface orientation in ARPES measurements, the Laue back-reflection x-ray 
measurements were performed at the ALS using a Realtime Laue delay-line detector from 
Multiwire Laboratories (Ithaca, NY) and diffraction spot indexing was performed using the 
MWL NorthStar orientation package. The x-ray Laue measurements along with the analysis 
results are presented in the Supplementary Figure SI 1. 

Spin-resolved measurements. Spin-resolved ARPES measurements were performed 
on the SIS beamline at the Swiss Light Source (SLS) using the COPHEE spectrometer 
with two 40kV classical Mott detectors and a photon energy of 24 eV, which systematically 
measures all three components of the spin of the electron (P x , P y , and P z ) as a function of 



its energy and momentum 



25 



26] . Energy resolution was better than 60 meV and 3% of 



the surface BZ for the spin-resolved measurements. Samples were cleaved in situ at 20 K 
at chamber pressure less than 2xl0 -10 torr. Typical electron counts on the detector reach 
5 x 10 5 , which places an error bar of approximately ±0.01 for each point on our measured 
polarization curves. 

Sample growth. Single crystals of Pb1_2.Sn3.Te were grown by the standard growth 



method (described in Ref. [39J for the end compound PbTe). High purity (99.9999%) 
elements of lead, tin, and tellurium were mixed in the right ratio. Lead was etched by 
CH3COOH+H2O2 (4:1) to remove the surface oxide layer before mixing. The mixture 
was heated in a clean evacuated quartz tube to 924 °C where it was held for two days. 
Afterwards, it was cooled slowly, at a rate of 1.5 "C-h^ 1 in the vicinity of the melting point, 
from the high temperature zone towards room temperature. The nominal concentrations 
were estimated by the Pb/Sn mixture weight ratio before the growth, in which the two 
representative compositions presented in the paper were estimated to be x nom = 0.20 and 
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^nom = 0.36. The concentration values were re-examined on the single crystal samples after 
growth using high resolution x-ray diffraction measurements by tracking the change of the 
diffraction angle (the 2Theta angle) of the sharpest Bragg peak, from which we obtained 
x = 0.22 and x = 0.39 (see Supplementary Figure S12 for diffraction data). In the paper we 
note the chemical compositions as Pbo.sSn .2Te and Pbo.6Sno.4Te. 

First-principles calculation methods. We perform first-principles calculations to 
extract both the electronic structure and the spin texture of the SnTe surface states. Our 
first-principles calculations are within the framework of the density functional theory (DFT) 
using full-potential projected augmented wave method [4Qj as implemented in the VASP 
package [41| . The generalized gradient approximation (GGA) |42[ is used to model exchange- 
correlation effects. The spin-orbital coupling (SOC) is included in the self-consistent cycles. 
The surfaces are modeled by periodically repeated slabs of 33-atomic-layer thickness, sep- 
arated by 13-angstrom-wide vacuum regions, using a 12 x 12 x 1 Monkhorst-Pack Appoint 
mesh over the Brillouin zone (BZ) with 208 eV cutoff energy. The room-temperature crystal 
structures of SnTe in ideal sodium chloride are used to construct the slab in order to have the 
required crystal structure and mirror symmetries required for the predicted TCI phase [ijj] 
(The rhombohedral distortion is found to only occur in the Sn-rich compositions (x > 0.5) 
and only at low temperature [l6]). The experimental lattice constants of SnTe with the 
value of 6.327 A are used j^j. The experimental lattice constant of PbTe with the value of 
6.46 A 43] in ideal sodium chloride lattice are used in the PbTe band structure calculation. 
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FIG. 1: Band inversion topological transition and surface states in Pbi_ z Sn x Te. a, 

The lattice of Pbi-^Sn^Te system is based on the sodium chloride crystal structure, b, The first 
Brillouin zone (BZ) of Pbi_ x Sn x Te. The mirror planes are shown using green and light-brown 
colors. These mirror planes project onto the (001) crystal surface as the X — T — X mirror lines 
(shown by red solid lines), c, ARPES measured core level spectra (incident photon energy 75 eV) 
of two representative compositions, namely Pb0.sSn0.2Te and Pbo.6Sno.4Te. The photoemission 
core levels of tellurium Ad, tin Ad, and lead 5d orbitals are observed, d, The bulk band-gap of 



Pbi_ x Sn x Te alloy system undergoes a band inversion upon changing the Pb/Sn ratio 8Hl0l|. A 



TCI phase with metallic surface states is theoretically predicted when the band-gap is inverted 
(toward SnTe) [l|]]. The lowest-lying conduction and valence states associated with odd and even 
parity eigenvalues are labeled by Lq and Lg , respectively. e,f, First-principles based calculation 
of band dispersion (e) and iso-energetic contour with energy set at 0.02 eV below the Dirac 
node energy (f) of the inverted end compound SnTe as a qualitative reference for the ARPES 
experiments. The surface states are shown by the red lines whereas the bulk band projections are 
represented by the green shaded area in e. 
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FIG. 2: Comparison of band insulator and topological crystalline insulator phases, a, b, 

ARPES low energy electronic structure measurements on Pb0.sSn0.2Te (x = 0.2). c, First-principle 
calculated electronic structure of PbTe (x = 0) as a qualitative guide to the ARPES measurements 
on Pbo.8Sno.2Te. PbTe, which lies on the non-inverted regime, is theoretically expected to be a 
conventional band-insulator [15]. d, e, ARPES low energy electronic structure measurements on 
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FIG. 2: (Previous page.) Pbo.6Sno.4Te (x = 0.4) under identical experimental conditions and 
setups as in (a and b). f, First-principle calculated electronic structure of SnTe (x = 1) as a 
qualitative guide to the ARPES measurements on Pbo.6Sno.4Te. SnTe, which lies on the inverted 
regime, is theoretically expected to be a topological crystalline insulater [la ]. The insets of panels 
a and d show the non-inverted and inverted bulk band-gap of the x = 0.2 and x = 0.4 samples, 
respectively. The blue and red lines in panels c and f are the calculated dispersion of the bulk 
bands and the surface states, respectively. 
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FIG. 3: Dirac surface states and topological crystalline insulator phase in Pbo.6Sno.4Te. 

a, First-principle calculated surface states of SnTe with energy set at 0.02 eV below the Dirac 
node energy are shown in red. The blue dotted lines mark the first surface BZ. b, Left panel: 
ARPES iso-energetic contour mapping [E-q = 0.02 eV) of Pbo.6Sno.4Te covering the first surface 
Brillouin zone (BZ) using incident photon energy of 24 eV (throughout Fig. 3). Right panel: 
Spectral intensity distribution as a function of momentum along the horizontal mirror line (defined 
by k y = 0). c, High resolution Fermi surface mapping (E-q = 0.0 eV) in the vicinity of one of the 
X points (indicated by the white square in b). d-f, Dispersion maps {E-q vs k) and corresponding 
energy (momentum) distribution curves of the momentum space 
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FIG. 3: (Previous page.) cuts 1, 2, and 3. The momentum space cut-directions of cuts 1, 2, and 
3 are denned by blue dotted lines in panel c. The second derivative image (SDI) of the measured 
dispersion is additionally shown for d. 
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FIG. 4: Spin polarization of topological Dirac surface states, a, ARPES dispersion map 
along the mirror line direction. The white dotted lines show the binding energies chosen for spin- 
resolved measurements, namely SR-Cut 1 at Eb = 0.06 eV and SR-Cut 2 Eb = 0.70 eV. Inset: 
Measured spin polarization profile is shown by the green and blue arrows on top of the ARPES 
iso-energetic contour at binding energy Eb = 0.06 eV for SR-Cutl. b,c, Measured in-plane spin- 
resolved intensity (b) and in-plane spin polarization (c) of the surface states (SR-Cut 1) near 
the Fermi level at Eb = 0.06 eV. d,e Measured in-plane spin-resolved intensity (d) and in-plane 
spin polarization (e) of the bulk valence bands (SR-Cut 2) at E-r = 0.70 eV. f, Out-of-plane spin 
polarization measurements of SR-Cut 1 and SR-Cut 2. The error bar is ±0.01 for data points in 
all spin polarization measurements, g, Theoretically expected spin polarization configuration of 
the surface states which corresponds to a mirror topological invariant 
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28j ] . Measured spin polarization 



FIG. 4: (Previous page.) (the mirror Chern number) um = —2 27, 
texture configuration (green and blue arrows) of the surface states (SR-Cut 1) is shown on top of 
the calculated surface states. The blue and red lines are the calculated bulk bands and surface 
states respectively. 
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FIG. 5: A comparison of Z2 (Kane-Mele) topological insulator and topological crys- 
talline insulator (TCI) phases, a, ARPES and calculation results of the surface states of a 
Z2 topological insulator GeBi2Te4, an analog to Bi2Se3 (4]. Left: ARPES measured Fermi surface 
with the chemical potential tuned near the surface Dirac point. Middle: First-principles calculated 
iso-energetic contour of the surface states near the Dirac point. Right: A stack of ARPES iso- 
energetic contours near the F point of the surface BZ. b, ARPES measurements on the Pbo.6Sno.4Te 
{x = 0.4) samples and band calculation results on the end compound SnTe. The green straight 
lines in both (a) and (b) represent the linear energy-momentum dispersion of the Dirac cones. 
The green circles in (b) represent the double Dirac cone contours near each X point on the surface 
of Pbo.6Sno.4Te. 
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Supplementary Figure 1: Calculation of bulk band dispersion in SnTe. Theoretical 
calculated electronic structure of SnTe is shown along r(first surface BZ) — X — f(second surface 
BZ) mirror line direction. The SnTe lattice in ideal sodium chloride structure without rhombohedral 
distortion is used in the calculation. The experimental lattice constants of SnTe with the value of 
6.327 A are used. Therefore the TCI surface states are realized in the calculation. The red and 
blue lines correspond to the TCI surface states and the bulk bands respectively. The dispersion of 
the surface states is found to be independent of the k z values, whereas the bulk bands are found 
to show strong k z dispersion with a k z period of An. Right panel: Theoretical calculation of SnTe 
surface states and projection of band bands, showing how the surface states are connected to the 
bulk bands. The bulk band projection is the bulk valence band integrated over the entire k z range 
from k z = to k z = An. 
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Supplementary Figure 2: Calculation of bulk band dispersion in PbTe. Theoretical 
calculated electronic structure of PbTe is shown along F(first surface BZ)— X — r(second surface 
BZ) mirror line direction. The experimental lattice constant of PbTe with the value of 6.46 A in 
ideal sodium chloride are used in the PbTe band structure calculation. The trivial insulator phase 
is realized in PbTe and no surface states is found in the calculation. The blue lines correspond 
to the bulk bands at different k z values. The bulk bands are found to show strong k z dispersion 
with a period of 4-7T. Right panel: Theoretical calculation of PbTe projection of band bands. The 
bulk band projection is the bulk valence band integrated over the entire k z range from k z = to 
k 7 = 4ir. 
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a pb o.s Sn o.2 Te (x=0.2) 
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6k = k - k(X) 

Supplementary Figure 3: Incident photon energy dependence study of Pb0.sSn0.2Te 
and Pbo.6Sno.4Te. a, The ARPES dispersion maps and corresponding momentum distribution 
curves (MDCs) along the f(first surface BZ) — X — f(second surface BZ) mirror line direction at 
different incident photon energies for the non-inverted composition Pb0.sSn0.2Te. b, The ARPES 
dispersion maps and corresponding momentum distribution curves (MDCs) along the T — X — T 
mirror line direction at different incident photon energies for the inverted composition Pbo.6Sno.4Te. 
For the non-inverted x = 0.2 sample, no surface states on the Fermi level is observed for all the 
incident photon energy values applied under the same experimental conditions and setups. The 
k z evolution (dispersion) of the bulk valence band at different k z values (different incident photon 
energies) is 
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Supplementary Figure 3: (Previous page.) in qualitative agreement with the theoretical calcu- 
lation results shown above. For the inverted x = 0.4 sample, a pair of metallic states is observed 
on the opposite sides of each X point. The dispersion of the metallic states on the Fermi level is 
found to show no observable change upon varying the incident photon energy The bulk valence 
band is observed to be a single hole-like band below the Fermi level at each photon energy, which 
shows strong dispersion with respect to the incident photon energy 
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Supplementary Figure 4: Spin measurements near the surface state Fermi level, a, 

ARPES measured dispersion map along the T — X — T mirror line direction. The white dotted 
line shows the binding energy chosen for the spin-resolved measurement of SR-Cut 1. SR-Cutl is 
at binding energy near the Fermi level (22b = 0.06 eV), and thus it measures the spin polarization 
of the surface states. Inset shows the momentum space direction of SR-Cut 1 in the (k x ,k y ) 
momentum space plane on top of the measured iso-energetic contour of the surface states at 
E-q = 0.06 eV. The blue and green arrows represent the measured spin polarization configuration 
of SR-Cutl. b-c, The spin-resolved intensity along ±y direction, and the y component of the 
measured spin polarization of SR-Cut 1. In total four spin vectors pointing along the ±y direction 
are clearly revealed by the spin-resolved measurements, d-e, and f-g, Same as b,c, but for the x 
and z components of the spin, respectively. No significant x and z components of spin polarization 
is observed. 
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Supplementary Figure 5: Spin measurements at high binding energy (bulk valence 
bands), a, ARPES measured dispersion map along the T — X — T mirror line direction. The white 
dotted line shows the binding energy chosen for the spin-resolved measurement of SR-Cut 2. SR- 
Cut 2 is at high binding energy far away from the Fermi level (E-q = 0.70 eV), and thus it measures 
the spin polarization of the bulk valence bands, b-c, d-e, and f-g The spin-resolved intensity 
and the measured spin polarization of SR-Cut 2 along the y, x, and z directions, respectively. No 
significant spin polarization is observed for SR-Cut 2 (the bulk valence bands). 
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Supplementary Figure 6: Surface state lineshape. a, Fermi surface mapping zoomed-in in 
the vicinity of an X point, b, Dispersion map, second derivative analysis, and the corresponding 
momentum distribution curves of the surface states along momentum space cut direction defined 
by the blue dotted lines in a. The momentum distribution curve at E-q = 0.25 eV is highlighted 
in red. c, The Lorentzian fitting of the momentum distribution curve at E-q = 0.25 eV. d, The 
dispersion data points k) of the surface states are obtained from the Lorentzian fitting of the 
momentum distribution curves at different binding enerigies and shown by the solid blue circles. 
Then the surface state dispersion is fitted by linear function in order to obtain the velocity of the 
surface states. 
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Supplementary Figure 7: Spin texture of the (001) SnTe surface states in the vicinity 
of an X point based on first-principles calculation. The in-plane components (P x ,P y ) 
of the SnTe surface states spin texture. The red arrows represent the spin polarization at the 
corresponding momentum space location. The direction and the length of the arrows show the 
direction and the magnitude of the spin polarization, respectively. The insets represent the length 
of the unit spin vector with spin polarization of 1 (100%). The out-of-plane component is found 
to be zero throughout the calculated (E-B,k x ,k y ) range. The two blue dots denote the two Dirac 
points on the opposite sides of the X point. The color scale shows the energy of the surface states' 
contours with respect to the Dirac point energy (Ed). The black dotted lines show the iso-energetic 
contours of the surface states at different energies. The left panel shows the spin texture below 
the Dirac points (lower Dirac cones), and the right panel shows the spin texture above the Dirac 
points (upper Dirac cones). 
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Supplementary Figure 8: Calculated spin texture decomposed into their orbital ori- 
gins. The spin of s, p x , p y , and p z orbitals of Sn and Te atoms are shown for the lower Dirac cones 
of the surface states. The inset represents the length of the unit spin vector with spin polarization 
of 1 (100%) for all panels. The two blue dots denote the two Dirac points on the opposite sides 
of the X point. The color scale shows the energy of the surface states' contours with respect to 
the Dirac point energy (Ed). The black dotted lines show the iso-energetic contours of the surface 
states at different energies. 
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Supplementary Figure 9: Orbital projections of the SnTe (001) surface states. The 

red lines denote that the corresponding electronic states have a strong contribution (weight) from 
certain orbitals, whereas the green lines denote that the corresponding electronic states have a very 
small contribution (weight) from certain orbitals. 



35 



Meas. Spin Structure 



b Crossing is suggested by 
two bands with opposite 
spins 




c Avoided crossing is 
possible for the same 
spins 




Momentum 



Momentum 



Mirror Line 




Supplementary Figure 10: Spin polarization profile and Dirac crossing, a, The spin 
polarization of the surface states along the T — X — T mirror line direction revealed by our spin- 
resolved measurement SR-Cut 1. b, Theoretically, two spin singlet bands with opposite spin 
directions can only cross each other, c, Avoided crossing is allowed for band with the same spin 
flavor, or a spin degenerate band. 
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Supplementary Figure 11: Scattering geometry and sample characterization, a, Top: 
The electron analyzer setup used to measure emitted photoelectrons from the sample surface. 
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Supplementary Figure 11: (Previous page.) The two hemispheres detect the electrons to pro- 
duce a 2D image of energy vs momentum in one shot. Bottom: The geometry of the detector 
relative to the sample surface. The momentum of the electron inside the sample can be extracted 
from the measured values of Ek; n , 9 and <p. b, X-ray Laue diffraction pattern of a representative 
Pbo.6Sno.4Te sample used for ARPES experiments. The Miller indices are noted for the diffrac- 
tion peaks in Laue pattern, which reveals that the cleavage surface of the crystal is perpendicular 
to the [001] vertical crystal axis, c, The 2Theta values of the (200) Bragg peak (the sharpest 



and most intensive peak for X-ray diffraction on Pbi_ z Sn x Te 



44 



451 ] . also see Supplementary 



Figure [T2k ) and the corresponding lattice constant values for PbTe (x = 0%), Pbp.74Snn 9fiTe 
(x = 26%), Pbo.50Sno.50Te (x = 50%), and SnTe (x = 100%) obtained from [4J-|4J] respectively. 
Linear fittings are applied to the data points, from which we obtain 2Theta = 0.00605a; + 27.594 
and a = -0.00133a; + 6.459. 
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Supplementary Figure 12: X-ray diffraction measurements on two representative com- 
positions. 
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Supplementary Figure 12: (Previous page.) a, X-ray diffraction measurements on the two 
representative compositions. The intensity of the XRD measurements is normalized (maximum 
intensity to be 1). The XRD data of Pbo.6Sno.4Te is offset by 1 with respect to that of Pb0.sSn0.2Te. 
The corresponding indices of the Bragg peaks are noted, b, The diffraction data zoomed-in in close 
vicinity of the (200) Bragg peak for the two compositions. The 2Theta values of the (200) Bragg 
peak are found to be 27.73 deg and 27.83 deg for the two compositions. Using the formula 2Theta = 
0.00605x + 27.594 obtained from linear fittings in Supplementary Figure [TTc. we obtain that these 
2Theta values correspond to x values of x = 0.22 (for the nominal concentration x nom = 0.20 
samples) and x = 0.39 (for the nominal concentration x nom = 0.36 samples). In the paper we note 
the chemical compositions as Pb0.sSn0.2Te and Pbo.6Sno.4Te. 
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Supplementary Discussion 

• The topological distinction between TI and TCI phases 

— TCI vs. TI / Even vs. Odd number of band inversions: It is important 
to note that the Pbi„ x Sn x Te system contains an even number of inversions [9]. 
Traditionally even number of band inversions are believed to be topologically 
trivial, and only odd number of inversions are believed to give rise to topological 
insulator states [2]. For this reason, although the Pbi-^Sn^Te system has long 
been known to contain an even number of inversions, they were believed to be 
trivial insulators, (see Ref. [2]). Since the much studied Z 2 TIs can only be 
realized via odd number of inversions, thus the even number inversions in the 
SnTe system determines that such system is topologically distinct from both the 
TI phase and the trivial band insulator phase, which means TCI is a new phase 
of topological matter. 

— Irrelevance of Time-Reversal Symmetry Protection: The topological dis- 
tinction between TCI and TI can be directly seen from our ARPES measured 
surface states. For example, Fig. 5 in the maintext shows an ARPES comparison 
between a Z 2 (Kane-Mele) TI GeBi 2 Te 4 and our new Pbo.6Sno.4Te samples. As 
shown in Fig. 5 of the maintext, GeBi 2 Te 4 belongs to the well-understood single 
Dirac cone topological insulator, which is topologically the same as the famous 
Bi 2 Se3 or Bi 2 Te3 systems. GeBi 2 Te4 has only one surface state (odd number), 
which is found to enclose the time-reversal invariant momenta (TRIM) or the 
Kramers point (the f point). For this reason, the single Dirac cone surface state 
in GeBi 2 Te 4 (also Bi 2 Se 3 , Bi 2 Te 3 , and all other Z 2 TI systems) is topologically 
protected by time-reversal symmetry. In sharp contrast, in the Pb .6Sn .4Te sys- 
tem, none of the surface states are found to enclose any of the time-reversal 
invariant momenta (TRIM). Thus the surface states in the Pb .6Sn .4Te system 
are irrelevant to time-reversal symmetry type of protection. The topological pro- 
tection in the case of Pbo.6Sno.4Te is a result of the mirror (spatial) symmetries 
of the crystal, as theoretically predicted by Ref. [15]. One relevant experimental 
evidence for the mirror protection is that all of the observed surface states are 
located on the mirror line momentum space directions (f — X — f ). 
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• Spin polarization, Dirac crossing and gapless nature 

Here we provide one experimental observation relevant to the Dirac crossing and gap- 
less nature: the spin polarization profile. As shown in Supplementary Figure [TOj for 
two electronic bands with opposite spin directions, crossing is the only fate. The 
avoided crossing and gap behavior is only allowed if the bands have the same spin 
flavor or they are spin degenerate. Here we show that at energy level 60 meV below 
the chemical potential (thus about 60 ± 20 meV around the Dirac point), the two 
branches of each surface state cone show the expected opposite spin directions. It has 
been recently experimentally shown that the spin texture of the Dirac surface states 
will be strongly disturbed when a gap is opened at the Dirac point [38]. The energy 
scale where the spin texture is disturbed roughly corresponds to the gap value [38]. 
Therefore, by showing a helical spin structure at about 60 ± 20 meV around the Dirac 
point in Supplementary Figure Hfc and Supplementary Figure [TOa . we can safely ex- 
clude a gap value of 50 meV. However, whether there is a < 10 meV gap or not in our 
surface states still need further studies. 
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Supplementary Methods 

• Laue measurements 

The surface cleavage termination is determined by X-ray diffraction measurements on 
the cleaved samples after ARPES measurements. Supplementary Figure [TTb shows a 
representative Laue measurement on the Pbo.6Sno.4Te (x = 0.4) sample. The Miller 
indices are noted for the diffraction peaks in Laue pattern, which reveals that the 
cleavage surface of the crystal is perpendicular to the [001] vertical crystal axis. The 
(001) surface termination determined here from Laue results is consistent with the 
four-fold symmetry of the ARPES Fermi surfaces and the size of the surface BZ in 
ARPES measurements. 

• X-ray diffraction measurements 

As shown in Supplementary Figure [12] for our XRD measurements on the two rep- 
resentative compositions, the 2Theta values of the (200) Bragg peak are found to be 
27.73 deg (for the nominal concentration x nora = 0.20 samples, noted as Pb0.sSn0.2Te 
in the maintext) and 27.83 deg (for the nominal concentration x nom = 0.36 samples, 
noted as Pbo.6Sno.4Te in the maintext). Using the linear fitting above, these 2Theta 
values correspond to x values of x — 0.22 (for the nominal concentration x nom = 0.20 
samples) and x = 0.39 (for the nominal concentration x nom = 0.36 samples). In the 
paper we note the chemical compositions as Pb0.sSn0.2Te and Pbo.6Sno.4Te. 

• Incident photon energy dependence study on Pbi_ a Sn a Te 

As a qualitative guide to the ARPES incident photon energy dependence measure- 
ments, here we show the first-principles calculated bulk bands k z dispersion of the 
two end compounds inverted SnTe (Supplementary Figure [1]) and non-inverted PbTe 
(Supplementary Figure [2]) respectively. 

Now we present comparative incident photon energy dependence studies between the 
inverted Pb .6Sn .4Te samples and the non-inverted Pb .sSn .2Te samples under the 
identical experimental conditions and setups (beamline, incident photon energy val- 
ues, incident light polarization, sample surface preparation procedure, etc.) as shown 
in Supplementary Figure [3j A pair of surface states without observable k z dispersion 
is observed in the inverted x = 0.4 samples but absent in the non-inverted x = 0.2 
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samples, no surface state on the Fermi level is observed for all the incident photon 
energy applied under the same experimental conditions and setups as shown in Supple- 
mentary Figure The k z evolution (dispersion) of the bulk valence band at different 
k z values (different incident photon energies) is in qualitative agreement with the the- 
oretical calculation results shown above. The incident photon energy range from 26 
eV to 10 eV shown in Supplementary Figure [3] corresponds to a wide k z range from 
1.87T to O.I71", which corresponds to going from the top of the bulk BZ to the center of 
the bulk BZ. 

• Fitting results of the experimental chemical potential with respect to the 
Dirac point energy 

We apply a linear fitting to the ARPES dispersion near the Fermi level along Cut 2 
in the maintext, which enables us to extract the velocity of the surface states, as well 
as the chemical potential (E-p) with respect to the Dirac point energy (Ed). 

— Supplementary Figure shows the Fermi surface mapping zoomed-in in the 
vicinity of an X point. The blue dotted line defines a momentum space cut 
direction (which is the same as Cut 2 in Fig. 3 of the maintext). 

— Supplementary Figure [6b shows the dispersion map and the corresponding mo- 
mentum distribution curves (MDCs) along the momentum space cut defined by 
the blue dotted line in Supplementary Figure [6K 

— The MDC is the ARPES measured intensity distribution along the momentum 
axis at a fixed binding energy E-q. For example, the white solid line in Supple- 
mentary Figure [6b is at binding energy E& = 0.25 eV and the corresponding 
MDC is highlighted in red in the MDC panel in Supplementary Figure [6b. 

— We fit the MDCs by two Lorentzian peaks. We take the MDC Lorentzian fitting 
at Eb = 0.25 eV as an example, as shown in Supplementary Figure [6b. The 
Lorentzian peak position reveals the momentum location of each branch of the 
Dirac surface states at E-q = 0.25 eV. 

— By performing such fitting at different binding energies, we obtain the dispersion 
of the surface states, namely (E B , k). The solid blue circles in Supplementary 
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Figure Oil show the obtained dispersion data points of the surface states from 
E B = 0.25 eV to E B = 0.05 eV with a binding energy step of AE B = 0.02 eV. 

— Then we fit the surface states' dispersion data points by a linear function. For 
the two branches of the surface state cone we obtain the following: for the branch 
with positive slope, we have E B = (2.68(±0.07)eV-A)-A; + 0.008(±0.009)eV. 
And for branch with negative slope, we have E B = (— 2.88(±0.08)eV-A)-A; — 
0.010(±0.008)eV. Thus, the linear fittings of both branches give the velocity 
(the absolute value of the slope) of the Dirac cone to be around 2.8 eV-A. The 
Dirac point energy Ed is given by the binding energy intercept (E B (k = 0)) of 
the linear fitting, which is 0.008 eV and -0.010 eV for each fitting respectively. 
This means the Dirac point energy Ed implied by the fitting of the positive slope 
branch is 0.008 eV below the chemical potential E-p, whereas the Dirac point 
energy Ed implied by the fitting of the negative slope branch is 0.010 eV above 
the chemical potential E-p. 

— In the maintext, we report a surface state velocity of 2.8 eV-A along this mo- 
mentum space cut direction (Cut 2 in Fig. 3 of the maintext). Furthermore the 
chemical potential with respect to the Dirac point is found to be E-p = Ed ±0.02 
eV 

• Spin-resolved measurements on Pbo.6Sno.4Te 

Here we show systematic spin dataset of these measurements for all three components 
of the spin polarization vector. The x, y directions for spin polarization vectors are 
defined in the inset of Supplementary Figure Hh- And the z direction follows the right- 
hand rule. The P x , P y and P z measurements for SR-Cut 1 and SR-Cut 2 are shown 
in Supplementary Figure H] and Supplementary Figure [51 As discussed in the main 
paper, in total four in-plane tangential (P y ) spins are revealed by these polarization 
measurements (Supplementary Figure Hfc) on the surface states. We note that the 
magnitude of the spin polarization is between 10% — 20%, which is smaller than the 
typically observed 40% net spin polarization of the Z 2 topological insulator surface 
states [481 ] . In fact the 10% — 20% spin polarization in our x = 0.4 samples is found to 
be consistent with the first-principles based spin texture calculation on SnTe surface 
states (see below). As shown in Supplementary Figure Ht, the polarization magnitude 
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of the inner two branches are found to be smaller as compared to that of the outer 
two branches. This can be understood by two independent reasons: First, since the 
inner two branches are close to each other in momentum space, their polarization 
signal can intermix with each other due to the finite momentum resolution of the 
spin-resolved measurements (the momentum resolution is about 3% of the surface BZ 
which corresponds to roughly 0.04A ). Second, the smaller polarization of the inner 
two branches is even observed in our calculation on SnTe (see below). This effect is 
related to the anisotropy of the surface states and the hybridization of the two adjacent 
Dirac cones. 

• Spin texture calculation of SnTe 

The spin texture of SnTe is calculated using VASP first-principles method. We have 
assumed the SnTe lattice in ideal sodium chloride structure without rhombohedral 
distortion. The calculated spin texture of the SnTe surface states is shown in Supple- 
mentary Figure [7] for the lower and upper Dirac cone, respectively. The out-of-plane 
component of the spin texture is found to be zero. We compare the calculated spin 
texture with the spin polarization measurements shown above: Along the F — X — T 
mirror line direction (k y = in Supplementary Figured), both theoretical calculations 
in Supplementary Figure [T^l (low Dirac cone) and spin measurements in Supplemen- 
tary Figure Hb reveal in total 4 spins with alternating directions in going along the 
mirror line. The magnitude of the spin polarization in calculation is found to range 
10% to 30%, which is consistent with the experiments in Supplementary Figure Hb. 
The inner two spins close to the X point are found to have smaller polarization as 
compared to the outer two, which is also in agreement with the experimental results 
in Supplementary Figure Hb. 

In order to better understand the spin texture, we also calculate the atomic orbital 
contribution of the spin texture. For the (001) surface, as shown in Supplementary 
Figure [91 



